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Abstract 

We investigate a relationship network of humans located in a metric space where 
relationships are drawn according to a distance-dependent probability density. The 
obtained spatial graph allows us to calculate the average separation of people in 
a very simple manner. The acquired results agree with the well-known presence of 
the small-world phenomenon in human relationships. They indicate that this feature 
can be understood merely as a consequence of the probability composition. We also 
examine how this phenomenon evolves with the development of human society. 

Key words: Small-world phenomenon, random networks, spatial graphs, 
convolution. 



1 Introduction 



In the 1960's, the american social psychologist Stanley Milgram examined 
how people know each other and introduced a quantity named the degree of 
separation D. It is the number of people needed to bind two chosen persons 
via a chain of acquaintances. E. g., if persons A and B do not know each other, 
but they have a common friend C, their degree of separation is D(A, B) = 1. 
Milgram measured the mean degree of separation between people in USA and 
found a surprisingly small value, (D) = 6. This gave another name to this 
phenomenon - "six degrees of separation". 

Nowadays, if a network has small average distance between its vertices together 
with a large value of the average clustering coefficient, we say that small 
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world phenomenon is present and such a network is called small world network 
(SWN). We often encounter SWP in random networks. 

Clearly "To be an acquaintance" is a somewhat vague statement. There are 
various possible definitions - e. g. shaking the other's hands, talking to each 
other for at least one hour, etc. Fortunately, results do not depend significantly 
on the specific choice, SWN was observed in all of those cases. However, the 
number six in the name of the phenomenon can not be taken literally. Actually 
it is just an expression for the number, which is very small compared to the 
size of the investigated population, which is taken to be 6 400 millions (the 
approximate number of people on the Earth) in this article. 

Nowadays, SWP is a well known feature of various natural and artificial ran- 
dom graphs [1]. Article citations, World Wide Web, neural networks and other 
examples exhibit this feature [2,3]. 

There are many ways to construct a SWN. Some models are rather mathemati- 
cal and do not examine the mechanism of the origin of a network. They impose 
some heuristic rules (e. g. [4,5,6]) instead. Other models look for the reasons 
for the introduced rules. This is much more satisfactory from the physicist's 
point of view. The first such model is known as "preferential linking" [7]. It 
is quite reasonable for cases like the growth of the WWW, where sites with 
many links to them are well known and in the future will presumably attract 
more links than poorly linked pages. 

In this work we focus on the random network of human relationships. It evolves 
in a very complicated manner, therefore it is very hard to impose some well 
accepted rules for its growth. Hence, we do not look for the time evolution 
of human acquaintances. Instead, we inquire a static case with the random 
network already developed. 

If the acquaintance between A and B is present, we link them with an edge. 
We obtain the random graph of human relationships in this way. We can in- 
troduce a metric to this network by assigning a fixed position in the plane to 
every person. In order to obtain analytical results, we assume a constant pop- 
ulation density. In particular, we suppose that people-vertices are distributed 
regularly and form a square lattice in the plane. With proper rescaling, the 
edges of unit squares in this lattice have length 1. Further we assume that 
the probability that two people know each other, depends on their distance 
by means of some distribution function. This model should keep some basic 
features of the real random network of human relationships. 
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2 The Mathematical Model 



Let's have an infinite square lattice where squares have sides equal to 1 and 
there is one person in every vertex. We label the probability that two peo- 
ple with distance d know each other Q(d). We assume homogeneity of the 
population, therefore this probability function is the same for every pair. 

Summation of Q(d) through all vertices leads to the average number of ac- 
quaintances for any person which we denote N\. Next, we assume that the 
function Q(d) is changing slowly on the scale of 1. Therefore, we can change 
summation to integration and obtain 

N A = r dx r dyQ(Jx* + yA. (1) 

J — oo J — oo v ' 

Our aim is to quantify the average degree of separation (D) for couples with 
the same geometrical distance equal to b. To achieve this, we choose two such 
people and label them A and B with positions ta = [0,0], fs = [b,0] (this 
particular choice will not affect our results substantially). 



3 An Analytical Solution 

Every person in the lattice can be located by its coordinates [x, y] . We will 
denote the distance between X and Y as dxY- Let us introduce a symbol ~ 
for the relation of acquaintances. This is a binary relation which is symmetric 
but not transitive. The probability that X knows Y is then P(X ~ Y) — 
Q(dxy) = Qxy- 

We name P(D) the probability that the degree of separation for A and B with 
distance b is equal to the number D.lf we want to find out the average degree 
of separation in our present network, (D), we need to know the probabilities 
P(D) for all different values of D. At the moment, only -P(O) is known, since 
apparently -P(O) = Q(dAB) = Q(b). 



A B 




Let's examine the degree of separation D = 2. This means that there are two 
other persons on the path between A and B. We denote their coordinates as 
r*i = (xi,yi) and r 2 = (x 2 ,y2)- For the presence of such a track, edges Al, 12 
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and 2B are needed together with edges A2, IB and AB missing (see picture 
above). Since their presence is independent, we have 



P(2) = J2 QA1Q12Q2B (l - Qa 2 ) (l - Qib) (l - Qab) ~ 

1,2 

~ If QaiQuQ2b(i - QA2)(l - Qib)(i - Qab) dndf 2 . (2) 
1,2 

where the summation runs through various placements of persons 1 and 2. 
The change of summation to integration is possible due to the fact that Q(d) 
is changing slowly on the scale of 1. 

Here we utilized the fact that in probabilities addition rule P(AUB) = P(A) + 
P{B) — ?(Afl B) we can neglect the last term since probabilities P(A), P(B) 
are small and P(A(lB) is of the higher order of smallness. Unfortunately due 
to this approximation we apparently reach "probabilities" P{D) higher than 
1 for high enough value of D. Though probabilities P(D) small with respect 
to 1 can be considered accurate. This implies that obtained results cannot 
be used to evaluate the exact value of the average degree of separation for 
nodes A and B because in such a calculation we would need value of P(D) 
for every D. Still from the growth of P(D) we can easily see for which D* 
it reaches relevant values, e. g. P(D*) = 1/3. This D* then characterizes the 
mean degree of separation of A and B. 

We can compute the first approximation to (2), getting 

P(2)<°> = JI QA1Q12Q2B dn df 2 . (3) 
1,2 

As Q A1 = (5(^i-0,yi-0), Q 12 = Q(x 2 -x 1 ,y 2 -yi), and Q 2B = Q(b-x 2 ,0-y 2 ) 
we notice that (2) is a double convolution of the function Q(d) enumerated at 
point (b, 0) . Thus we can write 

P(2)W=[Q*Q*Q](b,0) =}► P( J D) (0) =g [D] (6,0). 

For the Fourier transformation of the convolution, the following equation 
holds: 

^{Q [D] } = (nQ}) D - 

Using this formula we can write P(D)^ in the form 

P( jD )(°) = ^- 1 {(^[g]) D }(6,o). (4) 



The mean clustering coefficient (C) is the probability that two acquaintances 
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of A know each other. It can be evaluated in a way very similar to the calcu- 
lation of P(D), the corresponding graph is on the picture below. 




In order to write down an expression for (C) it is straightforward to rewrite 
(2). We obtain the number of connected triples A12 with node A fixed by this 
integration. We just have to avoid double counting of every track (interchang- 
ing positions of 1 and 2) - this brings an additional factor of 1/2. The average 
number of acquaintances for every vertex is Na, therefore the average number 
of possible triples is N A (N A -l)/2« N 2 J2. The mean clustering coefficient 
is the ratio of the average number of triples to the average number of possible 
triples. That is, 



(°) = ^2 // QA1Q12Q2A dfi df 2 = [Q * Q * Q] (0, 0) = 

A l2 A 



^^- 1 {(^[Q]) 3 }(o,o). 



(5) 



Equations (4) and (5) are solutions of the problem. Unfortunately, the relevant 
functions Q(d) (see next section) do not have an analytical form of their 
forward and inverse Fourier transformation. Therefore we have to calculate 
the values of (C) and P(D) numerically. Equation (4) requires a very high 
calculation precision. This makes the evaluation of P(D) very slow and even 
with some clever treatment (see Appendix A) it is in practise impossible for 
high values of b. This is just our case, because we are interested in b = 50 000. 
Thus some other (approximate) approach is needed. First we have to find more 
about the nature of function Q(d). 



4 An Empirical Entries 



In the present, there are approximately 6 400 millions people on the Earth. It 
means that length of the assumed square lattice side is 2L = 80 000. In order 
to obtain a numeric results we choose b = 50 000 and the average number of 
acquaintances iV A = 1 000. 

To get some insight on the distribution Q(d), some analysis is needed. First it 
is clear that Q(d) should be decreasing with d. Moreover, closely living people 
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know each other almost certainly. That is 



UmQ(d) = l. (6) 

a— >u 

Together with (1) we now have two requirements for Q(d). Indeed, there are 
many functions satisfying them, e. g. we can choose Q(d) = Cexp[— r/a]. 

The last quantity we can compute is the average number of distant people every 
person in the lattice know, N d . Here distant means that people's distance from 
the chosen fixed person (node) is greater than L/2. This is a simple analogy 
to the number of people we know on the other side of the Earth. So we have 

L/2 

N d = N A -2n J rQ{r) dr. (7) 
o 

If exponential distribution discussed above satistfies (1) and (6) it folllows that 
N d ~ 1CT 13 . This is in a clear contradiction to the fact that there are people 
who have very distant friends. Still we can improve N d if we use stretched 
exponential Q(d) = exp[— K d a ] with exponent a between 0.2 and 0.3. 1 How- 
ever, if we check Q(l) (probability to know our closest person) it is well below 
0.3. Stretched exponentials therefore satisfy condition (1) just formally and 
we will not it discuss it later. Moreover, mean clustering coefficient is then 
very small (from 2 • 10" 4 to 3 • 10~ 3 ). 

Now it's clear that distribution Q(d) can't decrease so fast as exponential 
functions, wide tails are inevitable in our model. This leads us to power-law 
distributions l/x a . According to (6) we demand 

QM = j^, (8) 

where b is fixed by (1). Number of far friends now ranges from N d ~ 0.01 
(a = 3.5) to N d ~ 17 (a = 2.5). This range of exponents gives us reasonable 
range for values of N d . 

In this article we also show results for the normal distribution Q n (d) = 
exp[— ad 2 } (N d ~ 0) and the uniform distribution within fixed radius Q u (d) = 
$(R A -d) (N d = 0). 

With regard to the fact that all used distributions Q(d) approach to zero for 
large values of d it is almost certain that the shortest chain of acquaintances 
between chosen A and B do not run out of the examined lattice with side 
80 000. Therefore it doesn't matter if we have integration (summation) bound 

1 Approximate solution presented in next section can be used also for this distribu- 
tion. 
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in infinity or ±L = ±40 000. This allows us to use all results derived for infinite 
lattice in the real case of finite lattice. 



5 An Approximate Solution for Power-law Distributions 



To demonstrate the calculation we take P(2) as an example again. In the previ- 
ous section we found out that power-law distributions are especially important 
in our model. Their joint probability Q{r\)Q(b — n) has sharp maximum for 
T\ — and low minimum for T\ = b/2. Their ratio is 

Q(b/2) 2 (Ay 



Q(b)Q(0) \b 

where a is the exponent in (8). This implies that in (2) we can constrain 
summation to Tai,Ta2 <C b or r B2 ,r B1 C 6 or r A i,r B2 <C b (see picture 
below). 



A B A B A B 

Here we obtained three different diagrams. Let's examine first one in detail. 
Since edges AB and Bl are long we can write 

P(2) » JJ Q A1 Q 12 Q 2B (l - Q A2 ) dfi df 2 . 

1,2 

It is easy to show that for power-law distributions Q(b — r-C) ~ Q(b) = Qab 
when r\ b. Thus we have (for corresponding diagram see picture below) 

P(2) ^ II Q A1 Q 12 Q AB (l - Q A2 ) dn dr 2 = 

1,2 

= Q(b) I f QaiQu dfi df 2 - Q(6) | f Q A1 Q 12 Q A2 dfi df 2 . 



1,2 1,2 



A P 

Both integrals are easy to compute. Second one brings average clustering co- 
efficient (C) into account. The result is 



P(2)*Q(b)N 2 A (l-(C)). 
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Remaining two diagrams for P(2) can be evaluated in the same way. 

In the computation of P(D) for higher values of D we encounter products 
of kind (1 — <3i3)(l — Q24) ■ ■ ■ even after neglecting probabilities Qij for long 
edges ij. Here we can make first order approximation 

(1-Ql 3 )(l-Q24)wl-Ql3-Q24 

which is valid almost everywhere except small spatial region that do not con- 
tributes substantially (see section Results and discussion). Moreover, second 
approximation considering terms Q13Q24 would increase evaluated probabili- 
ties. Therefore first approximation results will be some lower bound estimates 
oiP(D). 

Higher values of D introduce long closed loops of kind AY1 . . .nA (n < D). 
Corresponding integrals can be carried out in the same way like it was pre- 
sented in the derivation of (5). Finally we obtain 

Cn ^^II Q^Qv-Q"***" = ^ _1 {(^[^) n }(°' )- ( 9 ) 

A l2 A 

This helps us to find values of C n for any n. Clearly C2 = (C). With the use 
of such a closed loop integrals we can write 

P(0)=Q(6), 
P(l)=Q(b)N A 2, 

P(2)=Q(b)Nl(3 

P(3) = Q(b)Nl(A 

P(4) = Q(b)Ni(5 

P(5) = Q(b)N 5 A (6 



6 Results and Discussion 

In this section we summarize results for various distributions Q(d) ranging 
from the uniform Q u and normal Q n to power-law distributions Q3.5-Q2.5 (see 
(8)) and flat distribution Qer- This list is sorted according to the quantity of 
long shortcuts in such networks of relationships. 



-2C 2 ), 
-6C 2 -2C 3 ), 

- 12C 2 -6C 3 -2C 4 ), 

- 20C 2 - 12C 3 - 6C 4 - 2C 5 ) 
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Flat Distribution 



If we have flat distribution Qer, every pair of vertices is connected with the 
same probability p. It is shown in [4] that in the network consisting of N 
vertices holds 

m~ 111 N 

Here pN is the average number of acquaintances for a person in the network, 
pN = N A . We have N A = 1 000 and iV = 6 400 millions thus D^ R = (!)-l» 
2.3and(C)^0. 



Uniform Distribution Within Fixed Radius 



We can discuss such a case where every person knows just N A closest neigh- 
bors. This leads us to the distribution Q u (d) = $(R A — d) where distance R A 
is fixed by (1). It gives us -Ra = \J N A /n and therefore 

It's worth to note that we don't have any randomness in this model thus 



Normal Distribution 



The only distribution which allows us to evaluate (4) analytically is normal 
distribution Q Q . The result is 



P(D) 



D + l 



cxp 



Tib 2 



N A (D + 1) 



It was argued before that a solution of the equation P(D*) = 1/3 characterizes 
value of the mean degree of separation (D) n . For N A = 1000 and b = 50 000 
we can use some approximations which lead us to 



D* 



IT 



d: 



2 Randomness can be introduced by random placement of vertices. Hence we obtain 
so called random geometric graphs discussed in [8]. This approach is complementary 
to presented one where vertex placement is fixed but their connecting is due to some 
probability distribution. 
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The actual value of D* is about one third of D* (this is due to the existence of 
some longer connections in the network, although it is extremely suppressed 
by the exponential decay). We can note that both D* and D* scale with b 1 . 
This clearly differs from In b scaling of the Erdos-Renyi model. The clustering 
coefficient (C) can be evaluated easily both for normal and uniform distri- 
bution. We obtain high values of (C) (see graph below) in both cases. This 
agrees with our expectations. 



Power-law Distributions 



Numerical computation of coefficients C2, ■ ■ ■ , C5 with (9) is rather fast - their 
values are shown in the table below. 





a = 2.5 


a = 3.0 


a = 3.5 


c 2 


0.068 


0.154 


0.233 


c 3 


0.030 


0.090 


0.153 


c 4 


0.016 


0.059 


0.109 


c 5 


0.009 


0.042 


0.084 



Substituting these values into (10) leads us to values of mean degree of sep- 
aration which are marked in the Fig. 1. We see that power-law distributions 
Q(d) results in high values of mean clustering coefficients (C) = C 2 together 
with small values of D* (from 6 to 4). Thus small world phenomenon is clearly 
present in these networks. 

One can also ask for some comparison with the well-known Barabasi-Albert 
model. Mean vertex degree is then (k) = 2m and mean clustering coefficient 
is (C) — (m — 1) In 2 N/8N (here m is degree of just added vertices, see [9]). 
With respect to our choice N A = l 000, N = 6.4 • 10 9 it follows that m = 500 



InD* A 
8 

6 

4 

2 





(C) 
0.6 

0.4 

0.2 





n 3.5 3 2.5 ER 



H 1 1 1 +- 

n 3.5 3 2.5 ER 



Fig. 1. Graphs of mean degree of separation and clustering coefficients for various 
distribution functions. 
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and (C) 10 11 . This is certinaly nonrealistic value, our model gives better 
estimation of (C) . 

For presented values of coefficients Cj expressions in parentheses in (10) do 
not fall close to zero for quite wide range of values of D. Therefore we can 
(very approximately) write 

P(D) « Q{b)DN°. 

Solution of the equation P(D*) = 1/3 is approximately D* as a\nb/ \n Na- 
With b this scales as In b. This is very different from b l scaling of (D) for the 
uniform a normal distribution. Such a scaling is similar to the scaling in the 
Erdos-Renyi model, though values of clustering coefficient are kept high as we 
demanded in the introduction. 

Probability P(2) can be evaluated also by straightforward summation in ac- 
cordance with (2) although it takes huge amount of computer time. Obtained 
values agree very well with results presented above for all examined exponents 
but 2.5 - this case requires more computer time than it was given. Compu- 
tation of P(3) in the same way exceeds our computer possibilities for every 
exponent but we do not regard it necessary. 



Time Evolution and Some Limitations 



Human relationships in modern world are much more widespread than it was 
in the past. One can think of slowly changing exponent of the power-law 
distribution function Q(d) from large values to smaller (perhaps resulting to 
almost flat distribution in the future - internet helps to bridge the distances). 
According to the Fig. 1 we see that this would affect exact value of clustering 
coefficient. However it would remain high enough for wide range of exponents. 
Similarly changes of mean degree of separation are not important at all - it 
remains very small compared to the size of human population. 

Finally it has to be noted that in the described model we do not consider 
presence of some organized hierarchic structures in human society. E. g. chief 
of the firm knows his employees, but he also knows another chiefs who know 
their employees, etc. Amount of people involved in the hierarchical tree grows 
exponentially with the number of its levels. Such an arrangement therefore 
introduces additional way how to know each other with small resulting degree 
of separation. In presented calculation we didn't include this effect. Yet there 
is one important insight. If we proved the degree of separation being small 
without considering of the hierarchies, their presence would even decrease it. 
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7 Conclusion 



We have examined the mean degree of separation and the clustering coefficient 
for a random network of human relationships in this article. We were able to 
compute these quantities in our model. For a power-law decay of probability 
Q(d), we obtained a small mean degree of separation compared to the size of 
the network, along with a large value of the mean clustering coefficient. Both 
of these features are typical for small world networks. Thus we have shown 
that the small world phenomenon can be understood as a simple consequence 
of additivity of probabilities. 

We saw that the style of calculation depends on the used distribution Q(d). 
The computation was finished analytically for some special cases. In other 
cases, thanks to some approximations, we utilized the advantage of (5) where 
b do not enter the inverse Fourier transformation, making it easy to evaluate 
numerically. 

It's worth to note that the model solved herein is similar to the Watts and Stro- 
gatz model [6] where long shortcuts were introduced by a random rewiring pro- 
cedure. In our model long shortcuts are present thanks to wide tails of power- 
law distributions. This model model brings two basic advantages. First, the 
derivation and the resulting relations for C and D* are more simple. Moreover, 
our model has more realistic foundations. Nevertheless, the typical behavior 
of this model is the same as in previous models. The introduction of long 
shortcuts to the system decreases the average degree of separation rapidly, 
but also keeps the clustering coefficient high enough for the so called small 
world phenomenon to appear. 



A Numerics of the Fourier Transformation 

The Fourier integrals encountered in the solution of presented problem can 
not be solved analytically thus numerical techniques have to be used. In the 
inverse Fourier transform this is especially awkward because we meet rapidly 
oscillating term exp[ifru]. Here b is the distance between chosen persons A and 
B, by assumption big number (6 = 50 000). Therefore we have to compute 
Fourier transformation of f(d) very accurately. In order to make computation 
less demanding on the computer time, it is convenient to find some approxima- 
tion in the computing of the inverse Fourier transformation. We will continue 
with this derivation in the onedimensional case for the sake of simplicity. 

The Fourier transformation of the even function f(x) is an even real function. 
According to the (4) we are looking for the inverse Fourier transformation of 
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its n-th power, we will denote it g(u). It is also even real function. Therefore 
its inverse Fourier transformation is real function (sine-proportional terms 
vanish). Thus 



_^ oo ^ oo 

9(b) = 2^f g(u)e lbu du = — J g{u) cos[6w] du. 

—oo —oo 

This integral can be expressed as the sum of contributions from all periods of 
the cos[fru] function, I n = (2nn/b, 2n(n + (here n G N) 

oo ^ r 

9(b) = E S n( b ), s n(b) = — / g(u) cos[6w] du. 

n=-oo 27r / 

In the integrand of previous equation we can make Taylor expansion of g(u) 
around £ n = 2n(n + l/2)/6. Thereafter terms of kind -u m cos[6-u] emerge (m G 
N). Such integrals are easy to compute - first two terms of resulting expansion 
are then 

1 d 2 g 



S n (b) 



b 3 du 2 



it 2 — 6 d 4 (? 



Qb 5 du 4 



Finally we have 



9(b) 



If 5 



du 2 



7T 2 



Qb 5 



~ d^ 
du 4 



(A.l) 



This helps us to speed up inverse Fourier transformation - we do not have to 
know so many values of g(u). For every range I n evaluation of g(u) in three 
points (for numerical calculation of second derivative in the leading term of 
(A.l)) is enough. We just have to keep in mind that these points have to 
be close enough (with respect to 2n/b), otherwise we can obtain evidently 
incorrect results (e. g. g(b) = when border points have distance 2n/b). 



Acknowledgements 



The author would like to thank to staff of his department, especially to Martin 
Mojzis and Vladimir Cerny for valuable conversations and to Marian Klein for 
computer time. Acknowledgement belongs also to Jan Bod'a for introduction 
to the field, Miska Sonlajtnerova for her enthusiastic encouragement and my 
parents for their support. 



13 



References 



[1] Watts, D. J., Small worlds: The Dynamics of Networks between Order and 
Randomness, Princeton University Press (2003). 

[2] Dorogovtsev, S. N. and Mendes, J. F. F., Evolution of networks, Adv. Phys. 51 
(2002), 1079. 

[3] Reka A., Barabasi A.-L., Statistical mechanics of complex networks, Rev. of 
Modern Physics 74 (2002), 47. 

[4] Erdos P. and Renyi A., On random graphs, Publications Mathematicae 6 (1959), 
290. 

[5] Erdos P. and Renyi A., On the evolution of random graphs, Publ. Math. Inst. 
Gung. Acad. 5 (1960), 17. 

[6] Watts, D. J. and Strogatz, S., Collective dynamics of small-world networks, 
Nature 393 (1998), 440. 

[7] Barabasi, A. L. and Albert, R., Emergence of scaling in random networks, Science 
286 (1999), 509. 

[8] Penrose, M., Random Geometric Graphs, Oxford University Press (2003). 

[9] Fronczak, A., Fronczak, P., and Holyst, J. A., Mean-field theory for clustering 
coefficients in Barabasi- Albert networks, Phys. Rev. E 68 (2003), 046126. 



14 



